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We consider a qubit coupled to a nonlinear quantum oscillator, the latter coupled to an Ohmic 
bath, and investigate the qubit dynamics. This composed system can be mapped onto that of a qubit 
coupled to an effective bath. An approximate mapping procedure to determine the spectral density of 
the effective bath is given. Specifically, within a linear response approximation the effective spectral 
density is given by the knowledge of the linear susceptibility of the nonlinear quantum oscillator. To 
determine the actual form of the susceptibility, we consider its periodically driven counterpart, the 
problem of the quantum Duffing oscillator within linear response theory in the driving amplitude. 
Knowing the effective spectral density, the qubit dynamics is investigated. In particular, an analytic 
formula for the qubit's population difference is derived. Within the regime of validity of our theory, 
a very good agreement is found with predictions obtained from a Bloch-Redfield master equation 
approach applied to the composite qubit-nonlinear oscillator system. 

PACS numbers: 03.65.Yz,03.67.Lx,05.40.-a,85.25.-j 



I. INTRODUCTION 

The understanding of relaxation and dephasing prop- 
erties of qubits due to the surrounding environment is 
essential for quantum computation^. A famous model 
to study the environmental influences on the coherent 
dynamics of a qubit is the spin-boson model^Ti^, consist- 
ing of a two-level system (TLS) bilinearily coupled to a 
bath of harmonic oscillators. Although the bath degrees 
of freedom can be traced out exactly, analytical solu- 
tions are only possible vifithin perturbative schemes. First 
those perturbative in the coupling of the TLS to the bath 
are typically obtained within a Born-Markov treatment 
of the Liouville equation for the TLS density matrix^ 
or within the path integral formalisms,. The equivalence 
of both methods has been demonstrated restricting to 
low temperatures and low damping strengths in Ref. 0] . 
The second alternative approach is to perform perturba- 
tion theory in the tunneling amplitude of the two-level 
system. Within the so-termed non-interacting blip ap- 
proximation (NIBA)^^— it yields equations of motion for 
the TLS reduced density matrix enabling to capture the 
case of strong TLS-bath coupling. Reality is however of- 
ten more complex, as the qubit might be coupled to other 
quantum systems besides to a thermal bath. For exam- 
ple, to read-out its state, a qubit is usually coupled to a 
read-out device. 

In the following we mostly have in mind the flux qubit* 
read-out by a DC-SQUID. The latter mediates the dissi- 
pation originating from the surrounding electromagnetic 
bath and can be modeled both as a linear or nonlin- 
ear oscillator*"^^. Recently, the nonlinearity of qubit 
read-out devices, for example of a DC-SQUIDi^"— or 
a Josephson bifurcation amplifier—"—, has been used to 
improve the measurement scheme in terms of a faster 
read-out and higher fidelity. Specifically, the device was 
operated in a regime where the dynamics exhibited bifur- 



cation features typical of a classical nonlinear oscillator. 
As demonstrated e.g. in Ref. the quantum limit is 
within the experimental reach as well. 
From the theoretical side there are two different view- 
points to investigate the dynamics of a qubit coupled to 
an oscillator, with the latter in turn coupled to a thermal 
bath. The first way is to consider the TLS and the os- 
cillator as a single quantum system coupled to the bath, 
while the second is an effective bath description where 
the effective environment seen by the qubit includes the 
oscillator and the original thermal bath. The mapping to 
an effective bath has been discussed for the case in which 
the TLS is coupled to a harmonic oscillator in Ref. p3| . 
Specifically, the spectral density of the effective bath ac- 
quires a broadened peak centered around the frequency 
of the oscillator. This case has been investigated in Refs. 
[23 - l30j by applying standard numerical and analytical 
methods established for the spin-boson model. All those 
works showed that the peaked structure of the effective 
bath is essential when the eigenfrequency of the TLS be- 
comes comparable to the oscillator frequency. 
So far the first approach was used in Ref. [3l| to describe 
a qubit-nonlinear oscillator (NLO) system in the deep 
quantum regime. Here the effects of the (harmonic) ther- 
mal reservoir can be treated using standard Born-Markov 
perturbation theory. The price to be paid, however, is 
that the Hilbert space of the qubit-nonlinear oscillator 
system is infinite, which requires for practical calcula- 
tions its truncation invoking e.g. low temperatures^. 
In contrast to the above work we investigate here the case 
of a qubit-NLO system, with the latter being coupled to 
an Ohmic bath, within an effective bath description. Due 
to the nonlinearity of the oscillator, the mapping to a lin- 
ear effective bath is not exact. In this case a temperature 
and nonlinearity dependent effective spectral density well 
captures the NLO influence on the qubit dynamics. 
The paper is organized as follows: In section HIl we intro- 
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duce the model with the relevant quantities. In section 
mil the mapping procedure is investigated and the effec- 
tive spectral density for the corresponding linear case is 
given. Afterwards the mapping procedure is applied to 
the case of a qubit coupled to a nonlinear quantum oscil- 
lator. As a consequence of the mapping the determina- 
tion of the effective spectral density is directly related to 
the knowledge of the susceptibility of the oscillator. We 
show how the susceptibility can be obtained from the 
steady-state response of a quantum Duffing oscillator in 
section IIII CI In section IIVI the steady-state response 
of the dissipative quantum Duffing oscillator is reviewed 
and its susceptibility is put forward. The related effec- 
tive spectral density is derived in section |Vl In section |Vl] 
the qubit dynamics is investigated by applying the non- 
interacting blip approximation (NIBA) to the kernels of 
the generalized master equation which governs the dy- 
namics of the population difference of the qubit. A com- 
parison with the results of Ref. 31], obtained within the 
first approach, is shown. Further, analogies and differ- 
ences with respect to the linear case are discussed. In 
section IVIII conclusions are drawn. 



II. HAMILTONIAN 

We consider a composed system built of a qubit, -the 
system of interest-, coupled to a nonlinear quantum os- 
cillator (NLO), see Fig. [TJ To read-out the qubit state 
we couple the qubit linearly to the oscillator with the 
coupling constant g, such that via the intermediate NLO 
dissipation also enters the qubit dynamics. The Hamil- 
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FIG. 1: Schematic representation of the composed system 
built of a qubit, an intermediate nonlinear oscillator and an 
Ohmic bath. 



tonian of the composed system reads: 

Htot = Hs + i?NLO + -ffs+NLO + -ffNLO+B + Hb 

where 
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Here Hs represents the qubit Hamiltonian, where fj, is 
the particle's mass and U (q) a one-dimensional double 
well potential with minima at q = ±qo/2. -Wnlo is the 
NLO Hamiltonian, where the parameter a > accounts 
for the nonlinearity. When the oscillator represents a 
SQUID used to read-out the qubit, the oscillator fre- 
quency corresponds to the SQUID 's plasma frequency. 
The dissipation on the NLO is modeled in the following 
by coupling it to an Ohmic bath, characterized by the 
spectral density^: 



J(^) 



-6{lu- 



j) — rju! — MjLo. (3) 



In the classical limit it corresponds to a white noise 
source, where 77 is a friction coefficient with dimensions 
mass times frequency. 

In the following focus will be on the qubit dynamics in the 
presence of the dissipative nonlinear oscillator. Namely 
we will study the time evolution of the qubit 's position 
as described by: 



q{t) := Tr{ptot{t)q} ^Trs{prcd{t)q}, 



(4) 



where ptot and picd are the total and reduced density 
operators, respectively. The latter is defined as: 



Prcd 



TreTrNLolPtot}, 



(5) 



where the trace over the degrees of freedom of the bath 
and of the oscillator is taken. In Fig. [2] two different ap- 
proaches to determine the qubit dynamics are depicted. 
In the first approach, which is elaborated in Ref. (sij . 
one first determines the eigenstates and eigenvalues of 
the composed qubit-oscillator system and then includes 
environmental effects via standard Born-Markov pertur- 
bation theory. The second approach exploits an effective 
description for the environment surrounding the qubit 
based on a mapping procedure. This will be investigated 
in the next section. 



III. MAPPING 

The main aim is to evaluate the qubit 's evolution de- 
scribed by q{t). This can be achieved within an effec- 
tive description using a mapping procedure. Thereby 
the oscillator and the Ohmic bath are put together, as 
depicted in Figure [21 to form an effective bath. The ef- 
fective Hamiltonian 



(6) 



is chosen such that, after tracing out the bath degrees 
of freedom, the same dynamical equations for q{t) are 
obtained as from the original Hamiltonian iJtot- Due to 
the nonlinear character of the oscillator an exact mapping 
implies that HbqS represents a nonlinear environment. 
We shall show in the following subsection using linear 
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FIG. 2: Schematic representation of the complementary ap- 
proaches available to evaluate the qubit dynamics: In the first 
approach one determines the eigenvalues and eigenfunctions 
of the composite qubit plus oscillator system (yellow (light 
grey) box) and accounts afterwards for the harmonic bath 
characterized by the Ohmic spectral density J(ti;). In the ef- 
fective bath description one considers an environment built 
of the harmonic bath and the nonlinear oscillator (red (dark 
grey) box). In the harmonic approximation the effective bath 
is fully characterized by its effective spectral density Jcfi(cj). 



response theory that a linear approximation for HboS is 
justified for weak coupling g. Then Eq. © describes an 
effective spin-boson problem where 



Beff 
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and the associated spectral density is: 



A/" 



(8) 



The Hamiltonian ([6]) with ([7]) leads to coupled equations 
of motion^i^: 




TTijXj + rajiu^jXj — djq, 

where U'{q) = -^U{q). By formally integrating the sec- 
ond equation of motion and inserting the solution into 
the first equation the well-known Langevin equation for 
the operator q is obtained. This, in turn, allows to obtain 
the Langevin equation for qcs{t) :— Tr{/5cff(?(i)}^': 

^lqcs+^i^ dt'j^s{t-t')q^s + {U'{q))cS = 0, (9) 
Jo 

with the effective damping kernel 7off (i — t'). 
Notice that (. . . )cff indicates the expectation value taken 
with respect to pcff, which is the density operator asso- 
ciated to HeO^. In Laplace space, defined by 

yit) = I dXy{X)exp{Xt), (10) 



27ri 



yW 



dty{t) exp(— At), 



we obtain from Eq. ([9]) the equation of motion: 

^i\^q,s{\) + tiX-f,s{X)qcs{X) + (C/'(A))cff = 0. (11) 

The real part j'^ff{uj) = Re[7(A = —i^)] of the effective 
damping kernel 7cff(t) is related to the spectral density 



7off(w) 



Jeff('^) 



(12) 



The mapping for the case of zero nonlinearity a and 
Ohmic damping has been discussed in Ref. [2j|. There 
the influence of both the intermediate harmonic oscilla- 
tor and the bath is embedded into an effective peaked 
spectral density given by: 



tHO 



M(f]2 - a;2)2 + M-f'^uj^ 



(13) 



showing Ohmic low frequency behaviour J^{uj) 



Equation of motion for the nonlinear 
Hamiltonian 



As discussed above, the mapping requires the knowl- 
edge of the reduced dynamics of the system described by 
the variable q{t). Therefore we start from the coupled 
equations of motion derived from the Hamiltonian i?tot 
given in Eq. 



m + U'{q) = -gy, 
Mi) + riii + M VL^ii + aij^ 



(14a) 

-m + m. (14b) 



According to Eq. ([3]) , ry = Af 7 is the damping coefficient 
and 



m 



C0s{L0jt) + 



p 



(0) 



■ sin{ujjt) 



-Mjsit)m 

(15) 

a fluctuating force originating from coupling to the bath. 
In order to eliminate y from the first equation of motion, 
we have to calculate y[q{t)] from the second equation. 
In the following we look at equations of motion for the 
expectation values resulting from Eqs. (I14ap and (jl4bp . 
i.e., we look at the evolution of q{t) := Tr{ptot'?(t)} and 
y(t) := Tr{/5toty(t)}- Since we want to calculate y{t) 
we turn back to Eq. ([l} and treat the coupling term 
^^s+NLO as a perturbation, 'g ^ Mil^. Then the use of 
linear response theory in this perturbation is justified 
and we find: 



y{t) = (yW)o 



(16) 



+0{at) 



dt'e{t-t'){m,yit')])o9{q{t'))o0{t') 
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where (. . . )o denotes the expectation value in the 
absence of the coupling which we assume it has been 
switched on at time to =0. 

Notice that for a linear system, as for example the 
damped harmonic oscillator, the linear response becomes 
exact, such that the neglected corrections are at least 
of order 0{a'g^). Moreover, the time evolution of the 
expectation values is the same as in the classical case; 
this fact corresponds to the Ehrenfest theorem^. For 
nonlinear systems the expression in Eq. (1161) is an 
approximation, because all orders in the perturbation 
are nonvanishing^^ . 
In Laplace space, Eq. (fT6)) yields: 



(17) 



where Sy{X) = y{X) — {y{X))o and where xW is the 
Laplace transform of the response function or suscepti- 
bility: 



x{t-t') = --9{t-t'){[y{t),m])o- (18) 

Since q{X) - (g(A))o = 0{g^), from Eqs. I^T^ and (HH) 
it follows: 



tiX^q{X)+g\{X)q{X)+Oiafy) 
= -{U\X))-g{m)o- 



(19) 



That is, we have a normalization of the mass, and a 
damping-like term due to the coupled equations of mo- 
tion. The effect of the nonlinearity is embedded in the 
response function x- 

We assume in the following that in the absence of the 
coupling to the qubit the NLO and bath are in thermal 
equilibrium, which yields (y(i))o = for all times, and 
thus also: (y(A))o = 0. 



B. Mapping of the equations of motion and generic 
form for the effective spectral density 

By comparison of Eqs. pT|) and we can conclude 
that they yield the same dynamics if: 

(C/'(A))eff-(C/'(A)), (20) 

and the effective bath is chosen such that: 

^2X(A) 



/iA 



7cff(A). 



(21) 



By comparing the last equations with the relation p2 
and replacing A = — it follows: 



where x"('^) is the imaginary part of the susceptibility in 
Fourier space. We have now reduced the problem of find- 
ing the effective spectral density to that of determining 
the corresponding susceptibility. Notice that for a linear 
system the classical and quantum susceptibility coincide 
and are independent of the driving amplitude! In this 
case it is possible to calculate xi^) directly from the clas- 
sical equations of motion. For a generic nonlinear system, 
however, the classical and quantum susceptibilities differ. 



C. Linear susceptibility of a Duffing oscillator 

In order to evaluate the linear susceptibility, we solve 
the auxiliary problem of calculating the susceptibility of 
a quantum Duffing oscillator (DO), i.e., of the nonlin- 
ear quantum oscillator in Eq. ([2]) additionally driven by 
a periodic force with driving amplitude F and driving 
frequency ujex- The corresponding equation of motion is: 

M§ + T]ij + Mn^y + af = -Fe{t - ta) coa{uJext) + ^{t). 

(23) 

Application of linear response theory in the driving yields 
the equation for the expectation value of the position of 
the oscillator: 

y{t) = {m)o-i r dt'9{t^t'){[m,m])o 



xFcos{ujext') + 0{F^). 



(24) 



Jeff (a.) = -rx"i^), 



(22) 



Using the symmetry properties of the susceptibility x{^) 
we obtain in the steady-state limit: 

yst{t) = lim y{t) = {y{t))o + F cos{uJext)x'{^ex) 

to^ — oo 

+Fsm{i0ext)x"M+O{F^) 
= Acos{ujext + (j)) +0{F^). (25) 

Here the presence of the Ohmic bath implies 
limtp_j._oo(y(i))o = 0- Notice that due to symmetry in- 
version of the NLO, corrections of 0{F'^) vanish in Eq. 
(1^ . In Eq. A and (j) are the amplitude and phase of 
the steady-state response. It follows x(a;) = ^ exp(— i0), 
such that x"i^) — sin0. 



IV. STEADY-STATE DYNAMICS OF A 
DUFFING OSCILLATOR 

So far we have reduced the problem of finding the effec- 
tive spectral density to the one of determining the steady- 
state response of the Duffing oscillator in terms of the am- 
plitude A and the phase (j). These quantities were recently 
derived in Refs. [32l|33l|. using the framework of a Bloch- 
Redfield-Floquet description of the dynamics of the DO. 
The results in Ref. [33!] are applicable in a wide range 
of driving frequencies around the one-photon resonance 
regime ujex = + 3SyQ/(4?i) = fii for strong enough 
nonlinearities: y^F / [2^/2) <C 3ayg/4 ^ Ml, where we 
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introduced the oscillator length yo = \/ h/{MD,). 
As illustrated in Refs. [32ll33j the amplitude and phase 
are fully determined by the knowledge of the matrix el- 
ements of the stationary density matrix of the Dufhng 
oscillator in the Floquet basis, see e.g. Eqs. (67)-(70) 
in Ref. [s^. There the master equation yielding the 
elements of the stationary density matrix is analytically 
solved in the low temperature regime ksT << hfl impos- 
ing a partial secular approximation, yielding Eq. (70) of 
Ref. [33|, and restricting to spontaneous emission pro- 
cesses only. Here we follow the same line of reasoning as 
in Ref. [33[ to evaluate the amplitude and phase: we im- 
pose the same partial secular approximation and consider 
low temperatures fc^T < hfl. However, we include now 
both emission and absorption processes, i.e., we use the 
full dissipative transition rates as in Eq. (64) of Ref. [s^ . 
The imaginary part of the linear susceptibility x follows 
from the so obtained nonlinear susceptibility xn l in the 
limit of vanishing driving amplitudes: 



X'iuiex) = lim XNhi^^ex) 
F— >0 



(26) 



j/4J(J7i)2ni(0)4(2nth(ni) + 1)2 + - ^i) 



where 



ni(0) = 



1 



(27) 



For consistency also n\{{)) has to be treated up to first 
order in a only. 

Moreover, we used the spectral density J{lo) — M^lj and 
the Bose function nth{t) = [coth (Se/(2fcBr)) - 1] /2, 
which determines the weight of the emission and absorp- 
tion processes. 



V. EFFECTIVE SPECTRAL DENSITY FOR A 
NONLINEAR SYSTEM 
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FIG. 3: Comparison of the effective spectral density Jesiio) 
and a Lorentz curve for moderate damping. The param- 
eters are: Q, = 1.0, yoa/{hQ,) = 0.08, 7 = 0.097J7, and 
/3 = 10(fin)-\ 



VI. QUBIT DYNAMICS 

In the following we derive the dynamics of a qubit cou- 
pled to this effective nonlinear bath. Therefore we iden- 
tify the system Hamiltonian iJg introduced in Eq. ([T]) 
with the one of a qubit, denoted in the following as i?TLS- 
This is verified at low energies if the barrier height of the 
double well potential U{q) is larger than the energy sep- 
aration of the ground and first excited levels in each well. 
In this case the relevant Hilbert space can be restricted to 
the two-dimensional space spanned by the ground state 
vectors \L) and in the left and right potential well, 
respectively^. We start defining the actual form of the 
qubit Hamiltonian and its interaction with the nonlinear 
oscillator and afterwards introduce its dynamical quan- 
tity of interest, the population difference P{t). 



The effective spectral density follows from Eqs. (|22|) 
and ([M]) . It reads: 

Jcff(a;e^) (28) 

^ M-y^^ll{2nth{0,l) + l)2ni(0)4 + 4.MQ?{\u)^x\ - S^i)^ ' 

As in case of the effective spectral density , Eq. 
(|13p . we observe Ohmic behaviour at low frequency. 
In contrast to the linear case, the effective spectral 
density is peaked at the shifted frequency fii. Its shape 
approaches the Lorentzian one of the linear effective 
spectral density, but with peak at the shifted frequency, 
as shown in Fig. [H 

While in Refs. [32.33] the amplitude of the oscillator 
showed an antiresonant to resonant transition depending 
on the ratio of driving amplitude F and damping 7, the 
effective spectral density, obtained in the limit F ^ 
displays only resonant behaviour. 



A. Qubit 

The Hamiltonian of the TLS (qubit), given in the lo- 
calized basis {\L), is: 

^TLS = -^(£(7, + Act,), (29) 

where the corresponding Pauli matrices, 

the energy bias e accounts for an asymmetry between the 
two wells and A is the tunneling amplitude. The bias e 
can be tuned for a superconducting fiux qubit by applica- 
tion of an external flux ^cxt and vanishes at the so-called 
degeneracy point^i. For e A the states \L) and \R) 
are eigenstates of -ffxLS j corresponding to clockwise and 
counterclockwise currents, respectively. 
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The interaction in Eq. ([2]) is conveniently rewritten as 

--qoCTzyoia + a^) 



^^TLs-NLo = gqy (30) 

9 



:= hgaz{a + a^). 



Likewise we express the nonUnear oscillator Hamiltonian 
as: 



(31) 



= hn[j + -j+-^{a + a^) 



B. Population difference 



The dynamics of a qubit is usually characterized in 
terms of the population difference P{t) between the \R) 
and |L) states of the qubit: 



Pit) 



TrTLs{CTzPrcd(i)} 
(R|prcd(t)|R) - (L|/5red(i)|L), 



(32) 
(33) 



where /5rod(i) is the reduced density matrix of the TLS, 

Prcdit) = TTB{Pctf{t)}. (34) 

It is found after tracing out the degrees of freedom of 
the effective bath from the total density matrix Pcs{t) — 
exp^'s^"" * /5off (0) exp's:^'=ff*. it follows that in the two 
level approximation qcf[{t) — ^P{t), where qcs{t) is the 
position operator expectation value introduced in Sec. 

imi 

As we mapped the nonlinear system onto an effec- 
tive spin-boson model, the evaluation of the popu- 
lation difference P{t) of the TLS is possible using 
standard approximations developed for the spin-boson 
mode l^^i^^'^^ . Assuming a factorized initial condition 
Pcs{0) = /5TLs(0)exp(-^i?Boff/-^), the population dif- 
ference P{t) fulfills the generalized master equation 



Pit) 



dt'[K%t-t')P{t') + K''{t% t>0 (35) 



where K^{t — t') and K'^it — t') are symmetric and an- 
tisymmetric with respect to the bias, respectively. They 
are represented as a series in the tunneling amplitude. As 
an exact solution neither analytically nor numerically is 
available, due to the complicated form of the exact kernel, 
we impose in the following the so-called Non-Interacting 
Blip Approximation (NIBA)^i^. Applying NIBA corre- 
sponds to truncating the exact kernels to first order in 



and is therefore perturbative in the tunneling ampli- 
tude of the qubit. It is justified in various regimes: it is 
exact at zero damping, otherwise it is only an approx- 
imation which works at best for zero bias and/or large 
damping and/or high temperature^. One finds within 
the NIBA 



K'{t) = A'^cxp{-S{t))cos{R{t)), 
K^it) = A^exp{-S{t))sm{R{t)), 



(36) 



where S{t) and R{t) are the real and imaginary part of 
the bath correlation function: 



Q(t) = S{T)+lRiT)= [ 

Jo 



' Gcff(a;) 
duj ^—^ X 



(37) 



coth 



(1 — cos(a;t)) + i sm{ujt) 



where Gcs{oj) = q^Jcji{uj) / {T:h). In particu- 
lar, upon introducing the dimensionless constant 
= g'^^ni{Q)^ / {nVt^), we obtain: 



f + (|..|~f7i)2' 



(38) 



where we used riini(O)^ = f2 -|- 0(q;^), and 
7(;j := {2nth{^i) + 1)7/2. Consequently, the dynamics 
of the qubit is fully determined by the knowledge of 
the correlation function Q{t) and hence of the effective 
spectral density derived in section Ivl 
We now consider the qubit dynamics for the case of the 
effective nonlinear bath. Therefore we determine the 
actual form of the correlation functions S{t) and R{t). 
From Eq. ^ it follows: 



S{t) = Xt + L [exp(-7t;jT) cos (flir) - 1] 

+Zcxp(-7,;,r)sin(f]ir), (39) 
R{t) = I-exp{-j,^T)[Nsm{niT) (40) 

+IC0S {fiiT)] , 



where 



N 



Kith ^1 J 



(41) 



^^^^ 
/ 

L = - 



1 



Z = - 



1th cosh {f5hQ.i) - cos (^/i7th) 
[rJi sinh {l3hVLi) - 7t^ sin {Pfvyth)] 
I 1 



1th cosh {f5hQ.i) - cos ipn-ith) 
^th sinh {pmii) + VLi sin {Pfvyth)] ■ 

Here we have neglected the contribution coming from the 
Matsubara term, which is verified if the temperature is 
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high enough^, i.e. fcsT » frj/{2'7T). Moreover, we ap- 
phed in the contributions of the poles lying in the vicinity 
of ±f2i the approximation: 2J7i/(2r2i ± vJth) ~ 1- This 
corresponds effectively to neglect certain 0{^^f^) contri- 
butions. 



C. Analytical solution for the nonlinear peaked 
spectral density 

In this section we derive an analytical formula for the 
population difference P{t) for the symmetric case (e ~ 0), 
requiring weak damping strengths 7, such that a weak 
damping approximation of the NIBA kernels is verified, 
specifically 7/(27ril) << 1. As this calculation is ana- 
logue to the one illustrated in detail in Ref. [1^ , we only 
define the relevant quantities and give the main results. 
Due to the convolutive form of Eq. ([35| this integro- 
differential equation is solved by applying Laplace trans- 
form. In Laplace space it reads: 



P(A) 



(42) 



where P(A) = dtexp{—Xt)P{t) and analogously for 
K'^^'iX). 

Consequently, the dynamics of P{t) is determined if the 
poles of 



X + K'{X) = 



(43) 



are found and the corresponding back transformation is 
applied. We transform the kernels in Eq. ([55)) in Laplace 
space and expand them up to first order in the damping. 
This procedure is called weak damping approximation 
(WDA) in Ref. One obtains: 



K^''\X) 



K^^^X) 



dr exp(— At) exp(— S'o(t)) 



(44) 



{cos(i?o(r))[l - 5i(t)] - sin(i?o(r))i?i(T)} , 
0, 



where the indices {0, 1} denote the actual order in the 
damping. Specifically, 



Sir) - 5o(t)+5i(t)+0(72), 

R{t) = i?o(T)+i?l(T)+0(72), 



where 



So{r) = Y[cos{n,T) ^ 1], 

5*1 (r) = At cos{niT) + Bt + Csm{niT), 

i?o(T) = VFsin(r!ir), 

i?i(T) ^ V (l~ cos(riiT) - ^ sin(f7ir) 



(45) 
(46) 



(47) 



The zeroth order coefficients in the damping are given 
by: 



Y 
W 



-W 



1 



smh{l3hni) 
coshiPhfli] 
4,g^ni(0)4 

and the first order coefhcients by: 



(48) 



V 



Ph^i +sinh(^?il7i) 
cosliifihVLi) - I ' 
2.g2ni(0)S 



With this we are able to solve the pole equation for P(i), 
Eq. (|43p. as an expansion up to first order in the damp- 
ing around the solutions Xp of the non-interacting pole 
equation, i.e.. A* — Xp — "fKp + i"fv + 0{'^'^), as 7/O << I. 
Following Nesi et al<^ the kernel is rewritten in the com- 
pact form: 

/.oo 

X('')(A) = / drexp(-AT) {A^ ^cos(nf]ii) 



where 



[1 - ^i(t)] + A2 sin(nf}ii)i?i(T)} , (49) 



A„,c = Aexp(y/2)W(2-5„,o)(-j)"Jn(«o)cosh n 



A„ 



= Aexp(F/2)y'(2-5„,o)(-i)"J„(«o)sinh (^^^ 



and 



Mo 



Vni(0)4 



(50) 



(51) 



I 



To obtain analytical expressions, we observe that in the 
considered parameter regime where g/^l <^ 1 and j3Mli > 
1 it holds |uo| < 1. Following^* this allows effectively 
a truncation to the rt = and n = 1 contributions in 
K''-^\X) as the argument of the Bessel functions is small, 
leading to the following approximations: 



A^ - 



A^ 



A^ exp(r) Jo(uo) ~ A^ exp(r) (52) 
Vni(0)6- 



A^ exp 



02 



A^ exp(r) Vy^~~W^cosh(/3nrJi/2), 
4g2ni(0)6 



A^ 
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8 



Solving the undamped pole equation yields: 



(53) 



± 




A?,c 



IS? - 



Of 



-ni 



The last two equations allow to determine the oscillation 
frequency. Finally, within the WDA the qubit's popula- 
tion difference is obtained as: 



P{t) = exp(-7K-f) ^3 



-|-exp(-7K+f) 



A2. - \l_ 



cos 

cos(!:2+f) - sin(r2+t) 

i i-L 



(54) 



where k± = k(A±), which is derived in detail in Eq. 
(B.l.) of Ref. Note that for a consistent treatment 

if 0{g^) is kept we implicitly require 7 ^ <?, as only the 
first order in the damping is taken into account. 

We consider two possible resonance cases: First we 
choose the resonance condition fii = Aq^c, such that the 
oscillation frequencies are, to lowest order in Ai_c, 



= r^i =F 



(55) 



f^ + |«T5(l-|:«). 



As a consequence we obtain the so-called Bloch-Siegert 
shift: 



-n+ = 2g{l - 5a/2h), 



(56) 



which is also obtained in Ref. (31|. For comparison 
with^ we choose as second condition A = O, such that: 



3a Sag 
3a 3ag 



(57) 
(58) 



which agrees with the results of'^^ up to first order in the 
nonlinearity or /and in the coupling. Comparing witb^i 
we do not observe an exact agreement for the prefactors 
of the mixed terms of order 0{ag). 

We show in Figs. S] and [S] a comparison of the ana- 
lytic WDA formula Eq. ([M)) . the numerical solution 
of the NIBA Eq. ((351), denoted by NIB A, and the re- 
sults obtained in Ref. [Slj from a numerical solution 
of the Bloch-Redfield equations referred to as TLS-NLO 
approach. We observe that the dynamics is dominated 
by two frequencies and well reproduced within all three 
approaches. In the Fourier spectrum we observe tiny 
deviations of the resonance frequencies. There are two 
different reasons for these deviations: First the coupling 



strength g is large enough, that higher orders in the cou- 
pling yield a finite contribution in the effective bath de- 
scription. Second Eq. ([55| has to be expanded in both 
the nonlinearity and the coupling g, which is not possible 
in the numerical program. However, as we derived above 
when expanding the analytic formula, we find up to low- 
est order in the coupling g and in the nonlinearity a the 
same results. We emphasize that this small discrepancy is 
also seen for the corresponding linear system in the work 
of Hausinger et ali^ when comparing the NIBA results 
in Ref. [2§\ with those of the Bloch-Redfield procedure. 



p. 0.0 




50 100 

t[n-'] 



FIG. 4: Comparison of the behaviour of P{t) as obtained 
from the numerical solution of the Bloch-Redfield equations 
based on the TLS-NLO approach, Ref. [3l|, the numerical 
solution of the NIBA equation, Eq. (|35p . and the analytical 
formula provided in Eq. (|54|) . The chosen parameters are: 
a = 0.02(fin), g = O.lSfi, e = 0, 7/(27rn) = 0.0154 and P = 
10{hfl)-^. The dynamics agree within all three approaches. 




FIG. 5: Corresponding Fourier transform of P{t) as shown in 
Fig. SI 

To clarify the above statements we consider also the case 
that the coupling is weak, i.e., g <^ In the regime 

where the coupling is much weaker than the nonlinearity 
{Hg <C a), Eq. (32) given in Ref. {HI has to be expanded 
differently. Note that in this regime the results of Eq. 
(41) in Ref. [3l[ are not applicable. A proper expansion 
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allows in this regime to neglect 0{g'^) or higher if 0{a^) 
is neglected. The transition frequencies, when choosing 
f2 = A, are then determined by Eq. (32) of Ref. [31]: 



n, 



(59) 



Applying also an expansion of Eq. (|53p consistent with 
this parameter regime, we obtain: 



such that: 



1 



(60) 



(61) 



= ni = n + 3a/h. 



The transition frequencies in Eqs. ([59|) and (|6T|) coincide, 
and in Figs. [S] and [7] there is no deviation observed when 
comparing the three different approaches. 




t[n-'] 



FIG. 6: As in Fig. g] but for a smaller TLS-NLO coupling 
constant q = O.OOlSri. 
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FIG. 7: Corresponding Fourier transform of P{t) shown in 
Fig. El 



D. Influence on the qubit dynamics due to the 
nonhnearity- A comparison of the NIBA for hnear 
and nonlinear eff'ective spectral densities 

In this last section we want to address the effects of 
the nonhnearity onto the qubit dynamics. The compar- 
ison of linear versus nonlinear case is done at level of 
the numerical solution of the NIBA equation and shown 
in Figs. [8] and HI As already obtained in Ref. [sij . 
we observe that the transition frequencies are shifted to 
higher values compared to the linear case. As a conse- 
quence also the amplitudes associated to the transitions 
are modified. Moreover, we observe a decrease of the 
vacuum Rabi splitting compared to the linear case. Con- 
sequently, the effect of the nonhnearity of the read-out 
device can be observed in the qubit dynamics. 



B, 0.0 




50 100 

t[n-'] 



FIG. 8: P{t) within the NIBA when using the linear and 
the nonlinear effective spectral densities, Eas. l|13p and (I28p 
respectively. Parameters are: a = 0.02{hQ,) or a = 
respectively, g = 0.18^, £ = 0, j/{2-kQ,) = 0.0154 and 




FIG. 9: Corresponding Fourier transform of P{t) shown in 
Fig. [3 The effect of the nonhnearity is to increase the reso- 
nance frequencies with respect to the linear case. As a conse- 
quence the relative peak heights change. 
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VII. CONCLUSIONS 

In this work we determined the dynamics of a qubit 
coupled via a nonhnear oscillator (NLO) to an Ohmic 
bath within an effective bath description. We investi- 
gated an approximate mapping procedure based on lin- 
ear response theory, which is applicable in case of weak 
nonlinearities and small to moderate qubit-NLO cou- 
pling. We determined the effective spectral density in 
terms of the qubit-oscillator coupling and the linear sus- 
ceptibility of a nonlinear oscillator. The susceptibility 
was calculated for practical purposes from the periodi- 
cally driven counterpart of the original nonlinear oscil- 
lator. The so obtained spectral density shows resonant 
behaviour, specifically almost a Lorentzian form for the 
considered parameter regime, and is peaked at a shifted 
frequency, namely at the one-photon resonance between 
ground state and first excited state of the nonlinear oscil- 
lator. Moreover, this effective spectral density acquires 
a temperature dependence and behaves Ohmic at low 
frequencies. Based on the effective spectral density the 
qubit dynamics are investigated within the NIBA approx- 
imation. In addition an analytical formula for the qubit 
dynamics is provided, which describes very well the dy- 
namics at low damping. These results were compared 
to the numerical solution in Ref. [3l|, where the Bloch- 
Redfield equations for the density matrix of the coupled 
qubit nonlinear oscillator system (TLS-NLO) are solved. 



We find an overall agreement of the two approaches and 
show that deviations are of order 0{ag), where a is 
the nonlinearity and g the coupling strength. Exem- 
plarily this effect was analyzed for two possible coupling 
strengths g. We emphasize that parameters like tem- 
perature and damping and especially the strength of the 
coupling g and nonlinearity a determine the appropriate 
form of the expansions in the different regime of param- 
eters. Due to the tunability of the parameters various 
qubit dynamics are possible. In agreement with Ref. (3l| 
we observed the following effects due to the nonlinearity: 
First the transition frequencies of the two dominating 
peaks, where we are in the regime g ^ a/H, are shifted 
to larger values compared to the linear case. As a conse- 
quence also the amplitudes of the coherent oscillations of 
the population difference P{t) are modified. Moreover, 
the Bloch-Siegert shift is decreased due to the nonlinear- 
ity. 

We conclude that, as in case of the corresponding linear 
syste m^^''^^ , the effective bath description provides an al- 
ternative approach to investigate the complex dynamics 
of the qubit dissipative NLO system. 
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